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Low-frequency quadrupole vibrations in deformed 36,38,40 Mg are studied by means of the deformed 
Quasiparticle-RPA based on the coordinate-space Hartree-Fock-Bogoliubov formalism. Strongly 
collective K n = + and 2 + excitation modes (carrying 10-20 W.u.) are obtained at about 3 MeV. It 
is found that dynamical pairing effects play an essential role in generating these modes. It implies 
that the lowest K* = + excitation modes are particularly sensitive indicators of dynamical pairing 
correlations in deformed nuclei near the neutron drip line. 



in 
o 
o 



O 



> 
O 

o 
in 
o 



o 

=5 



X 



PACS numbers: 21.60.Ev, 21.60.Jz, 21.10.Re 



1. Introduction 

The physics of drip-line nuclei is one of the current 
frontiers in nuclear structure physics. The number of 
unstable nuclei experimentally accessible will remark- 
ably increase when the next generation of radioactive 
ion beam facilities start running, and we shall be able 
to study the properties not only of the ground states but 
also of low-lying excited states of drip-line nuclei in the 
medium-mass region. In order to quest for the new kinds 
of excitation mode unique to unstable nuclei associated 
with new features such as neutron skins, many attempts 
have been made using the self-consistent RPA based on 
the Skyrme-Hartree-Fock (SHF) method Q, Q and the 
Quasiparticle RPA (QRPA) including the pair correla- 
tions 0, 0, IS IS • Most of these calculations, however, 
are restricted to spherical nuclei, and low-frequency RPA 
modes in deformed unstable nuclei remain largely unex- 
plored, except for some recent attempts 0, H, [nj . 

InRef. EH, 

we investigated properties of octupole ex- 
citations built on superdeformed states in neutron-rich 
sulfur isotopes by means of the RPA based on the de- 
formed Woods-Saxon potential in the coordinate-mesh 
representation. We found that low-lying states created 
by excitation of a single neutron from a loosely bound 
low- SI state to a high- SI resonance state acquire extremely 
strong octupole transition strengths due to very extended 
spatial structure of particle- hole wave functions. We have 
extended this work to include pairing correlations self- 
consistently, and constructed a new computer code that 
carries out the deformed QRPA calculation on the basis 
of the coordinate-space Hartree-Fock-Bogoliubov (HFB) 
formalism |T^ . In this paper, we report a new result of 
the deformed QRPA calculation for low-frequency soft 
quadrupole modes with = 0+ and 2+ in 36,38,40 Mg 
close to the neutron dri p li ne. According to the Skyrme- 
HFB calculations 0, 0], these isotopes are well de- 
formed. The shell-model calculation [l{| also suggests 
that the ground state of 40 Mg is dominated by the neu- 
tron 2p-2h components indicating breaking of the N = 28 
shell closure. We have studied properties of soft modes 
of excitation in these Mg isotopes simultaneously taking 
into account the deformed mean-field effects, the pair- 



ing correlations, and excitations to the continuum, and 
made microscopic analysis of the mechanism of gener- 
ating these collective modes. The result of calculation 
suggests that the soft K * = + modes are particularly 
sensitive indicators of dynamical pairing correlations in 
deformed nuclei near the neutron drip line. 

2. Method 

2.1. Mean- field calculation 

Assuming axially symmetric deformation, we solve the 
HFB equation [H Ej, 



h"(ra) -(fc«(r£r)-A)Av2,o(» v ) 



■E r 



t Pl, a ( r(7 ) 

on a two-dimensional lattice discretizing the cylindrical 
coordinates (p, z). For the mean-field Hamiltonian h, we 
employ the deformed Woods-Saxon potential with the 
parameters used in [XjJ except for the isovector potential 
strength depending on the neutron excess, for which a 
slightly smaller value, 30 MeV, is adopted in order to 
describe 40 Mg as a drip- line nucleus in accordance with 
the Skyrme-HFB calculations pi Hi]. The pairing field h 
is treated self-consistently using the surface-type density- 
dependent contact interaction, 



Vpp (r,r') = V (l-?-^)6(r-r% 



(2) 



with Vn = -450 McV-fm 3 and 



g is (r) denotes the isoscalar density. We use the lattice 
mesh size Ap = Az = 0.8 fm and the box boundary 
condition at yO ma x = 10.0 fm and z maK — 12.8 fm. The 
quasiparticle energy is cut off at 50 MeV. 

2.2. Quasiparticle-RPA calculation 

Using the quasiparticle basis obtained in the previous 
subsection, we solve the QRPA equation in the standard 



0.16 fm 3 , where 
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matrix formulation, 
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(3) 

For residual interactions, we use the Skyrme-type in- 
teraction for the particle-hole channel, 

t. 



toil + xoP^ + ^il + xsP^g^ir) 8(r-r') 



(4) 

with t = -1100 MeV-fm 3 , t 3 = 16000 MeV-fm 6 ,a; = 
0.5, X3 — 1.0, P a being the spin exchange operator. 
The density-dependent contact interaction (J2J) is used 
for the particle-particle channel. Because the deformed 
Wood-Saxon potential is used for the mean-field, we 
renormalize the residual interaction in the particle-hole 
channel by multiplying a factor f p h to get the spurious 
= 1 + mode (associated with the rotation) at zero 
energy (v p h — * f P h ■ Vph)- We cut the model space for the 
QRPA calculation by E a + E (j ^ 30MeV which is smaller 
than that for the HFB calculation. Accordingly, we 
need another self-consistency factor f pp for the particle- 
particle channel. We determine this factor such that the 
spurious K * = + modes (associated with the number 
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fluctuation) appear at zero energy (v pp 

The intrinsic matrix elements (X\Q2k\0) of the 
quadrupole operator Qik between the excited state |A) 
and the ground state |0) are given by 



(A|Q a jf|Q) 
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where 
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dpdzQ^{p, z) 



We calculate the intrinsic strength functions, 

s 1s (lo) = J2 \(MQ¥k\Q)\ 2 s(^ - 



(6) 



(7) 



for the isoscalar quadrupole operators Q\ S ki an< i use the 
notation B(Q IS 2) = |(A|Q^|0)| 2 . 

3. Results and discussion 

Calculated quadrupole transition strengths for the 
K 7 * = 0+ and 2+ excitations in 36,38,40 Mg are displayed 
in Fig. n We see prominent peaks at 2.9 and 3.3 MeV 
associated with the K v = 2 + and + modes, respectively, 
in 36 Mg; at 2.8 MeV for the K 7 * = 2+ mode in 38 Mg; at 
2.8 and 2.9 MeV for the K 7 " = 2+ and 0+ modes, respec- 
tively, in 40 Mg. This fi gurc indicates that the transition 
strengths are very large (10-20 W.u.) and increase as we 
approach the neutron drip line. Figure |21 shows the ra- 
tios of the matrix elements for neutrons and protons. It 
is seen that the neutron excitation becomes dominant as 
the neutron number increases. 




12 3 4 
hoi (MeV) 

FIG. 1: Isoscalar quadrupole strengths B(Q IS 2) for the K = 
+ (upper panel) and K = 2 + excitations (lower panel) on the 
prolate ground states of 36,38,40 Mg, obtained by the deformed 
QRPA calculation with = 0.28. The arrows indicate the 
neutron threshold energy; Eth = 4.01MeV (lqp continuum), 



4.1 MeV (2qp continuum) for 
0.60 MeV for 40 Mg. 



3 Mg, 2.14 MeV for J8 Mg and 



Let us focus our attention on the properties of excited 
states in 40 Mg. Details of the transition strength for 
the IC = + excitations are displayed in Fig. [3] We 
immediately notice that the transition strength for the 
lowest excited state is significantly enhanced from that 
of the unperturbed two-quasiparticle excitations. From 
the QRPA amplitude listed in Table HI it is clear that 
this collective mode is generated by coherent superpo- 
sition of neutron excitations of both particle-hole and 
particle-particle types. Let us discuss the reason why 
this mode acquires eminently large transition strength. 
In Ref. |ll| . we pointed out several examples where a 
neutron excitation from a loosely bound hole state to 
a resonant particle state brings about very large transi- 
tion strength. This is a natural consequence of the fact 
that their wave functions are significantly extended from 
the nuclear surface. In the present calculation for 40 Mg, 
this effect is much enhanced due to the pairing correla- 
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FIG. 2: Ratios of the neutron and proton matrix elements 
Mv/Mk, divided by N/Z, are plotted as functions of neutron 
number. Filled circles and squares indicate the ratios for the 
lowest K w = + and 2 + modes, respectively. 
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TABLE I: QRPA amplitudes of the soft = 0+ mode at 2.9 MeV in 40 Mg. This mode has B(E2) = 2.4 e 2 fm 4 , B(Q"2) = 112 
fm 4 , and B(Q IS 2) = 148 fm 4 . The single-particle levels are labeled with the asymptotic quantum numbers [JVnsAJfi of the 
dominant components of the wave functions. Only components with \f a p\ > 0.1 are listed. 
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FIG. 3: Left: Isoscalar quadrupole strengths B(Q IS 2) for the 
K* — + excitations in 40 Mg are plotted in the upper panel, 
while unperturbed two-quasiparticle strengths are shown in 
the lower panel. The arrow indicates the threshold energy 
2 1 A | = 0.60 MeV. Right: Two-quasiparticle excitations gener- 
ating the soft K n — + mode at 2.9 MeV. The pertinent levels 
are labeled with the asymptotic quantum numbers [JVnaAJH. 
The Fermi surface for neutrons is indicated by the dashed 
line. 



tions among these loosely bound and resonance states. 
To show this point, we plot in Fig. ^spatial distributions 
of two-quasiparticle wave functions generating the soft 
K 77 = + mode. We see that these are significantly ex- 
tended beyond the half-density radius. In addition, their 
spatial structures are rather similar; this is a favorable 
situation to generates coherence among them. In Fig. [5] 
we show how the transition strengths for the IC = + 
and 2 + modes change if the residual particle-particle in- 
teractions are ignored. We see that, although the promi- 
nent peak for the K 71 = 2 + mode still remains, the tran- 
sition strength decreases. For the K 7 ' = + mode, we 
notice a more striking effect; the prominent peak seen 
in Fig. ^ disappears when the dynamical pairing corre- 
lation is absent. Therefore we can say that the coher- 
ent superposition of particle-hole, particle-particle and 
hole-hole excitations is indispensable for the emergence 
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FIG. 4: Spatial distribution functions Q^oiPt z ) f° r some two- 
quasiparticle excitations generating the soft K* — + mode 
in 40 Mg. The contour lines are plotted at intervals of 0.002. 
The thick solid line indicates the neutron half density radius. 



of this mode. The importance of the coupling between 
the (particle-hole type) f3 vibration and the (particle- 
particle and hole-hole type) pairing vibration has been 
well known in stable deformed nuclei [T3] • A new feature 
of the soft K n = + mode in neutron drip-line nuclei is 
that this coupling takes place among two-quasiparticle 
excitations that are loosely bound or resonances, so that 
the transition strengths are extremely enhanced. 

We examined the numerical stability of this soft K n = 
+ mode against variation of the box size. In Fig. 
we show the isoscalar and unperturbed two-quasiparticle 
transition strengths for the if" = + excitations in 40 Mg, 
obtained by using a larger box; p max x z max = 13.2 fm 
xl6.0 fm. We see that the position and the strength 
of the prominent peak are almost unchanged. This re- 
sult confirms that the main two-quasiparticle compo- 
nents generating this mode are weakly bound and reso- 
nant states, and that the numerical error associated with 
the continuum discretization is not very significant. 
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FIG. 5: Isoscalar quadrupole strengths for the K 71 = + and 
2 + excitations in 40 Mg, obtained by switching off the residual 
particle-particle interactions. 





180 n 


E 


150 




120 - 


JZ 

a> 


90 - 


c 


60 - 


CO 


30 - 




- 



isoscalar 



T 
1 



60 
50 ] 
40 
30 
20 
10 
^ 
2 3 4 5 
ha (MeV) 



unperturbed 



1 2 3 4 5 



FIG. 6: Left: The isoscalar transition strengths for the K n 
u excitations in 40 Mg calculated using the box with size 
Pmax x z max = 13.2 fm x 16.0 fm. Right: Unperturbed two- 
quasiparticle strengths. 



4. Conclusion 

We have investigated properties of excitation modes in 
deformed Mg isotopes close to the neutron drip line by 



means of the deformed QRPA based on the coordinate- 
space HFB. For 36,38,40^^ we i, ave obtained the soft 
= + and 2 + modes possessing large transition 
strengths. The microscopic mechanism of enhancement 
of the K v = + strength is especially interesting. There 
are two causes for this enhancement; 1) the strong cou- 
pling occurs between the quadrupole shape vibration and 
the neutron pairing vibration, and 2) the two quasipar- 
ticle configurations coherently generating the K* = + 
modes possess significantly extended spatial structures. 
The first cause has been well known in studies of sta- 
ble deformed nuclei p"?|. but the second cause is unique 
to unstable deformed nuclei near the neutron drip line. 
The nature and roles of pair correlations for collective 
excitations in drip-line nuclei is currently under lively 
discussions jig. On the basis of the deformed QRPA 
calculation demonstrating that the soft K* = + mode 
is quite sensitive to the dynamical pairing, we suggest 
that an appearance of this soft mode is a very good indi- 
cator of pairing correlations in deformed drip-line nuclei. 
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